Soft wall effects on interacting particles in billiards 
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The effect of physically realizable wall potentials (soft walls) on the dynamics of two interacting 
particles in a one-dimensional (ID) billiard is examined numerically. The ID walls are modelled by 
the error function and the transition from hard to soft walls can be analyzed continuously by varying 
the softness parameter a. For a — > the ID hard wall limit is obtained and the corresponding wall 
force on the particles is the 5-function. In this limit the interacting particle dynamics agrees with 
previous results obtained for the ID hard walls. We show that the two interacting particles in the 
ID soft walls model is equivalent to one particle inside a soft right triangular billiard. Very small 
values of a substantiously change the dynamics inside the billiard and the mean finite-time Lyapunov 
exponent decreases significantly as the consequence of regular islands which appear due to the low- 
energy double collisions (simultaneous particle-particle-ID wall collisions). The rise of regular islands 
and sticky trajectories induced by the ID wall softness is quantified by the number of occurrences 
of the most probable finite-time Lyapunov exponent. On the other hand, chaotic motion in the 
system appears due to the high-energy double collisions. In general we observe that the mean finite- 
time Lyapunov exponent decreases when o increases, but the number of occurrences of the most 
probable finite-time Lyapunov exponent increases, meaning that the phase-space dynamics tends 
to be more ergodic-like. Our results suggest that the transport efficiency of interacting particles 
and heat conduction in periodic structures modelled by billiards, will strongly be affected by the 
smoothness of physically realizable walls. 
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I. INTRODUCTION 

Although physically realizable potentials arc inher- 
ently soft, most billiard models used in the literature have 
hard walls. For example, the Sinai billiard [l[,the Buni- 
movich stadium Q or the Annular billiard Q, among 
others, have hard walls and were used successfully to 
study the fundamental properties of classical and quan- 
tum chaotic systems. In such models the chaotic motion 
of the single particle dynamics arises as the consequence 
of the spatial billiard geometry. The question now is 
about the effect and importance of physically realizable 
potentials on the particles dynamics inside the billiards. 
Some works in this direction have shown that introducing 
soft walls do not destroy trajectories found in the hard- 
wall limit Q and may induce the appearance of regular 
islands in phase space @, S 0] ■ Such regular islands in- 
side the chaotic sea induce a "sticky" (or trapped) motion, 
which is a common phenomcnum in conservative systems 
Q. They arise from broken Kolmogorov-Arnold-Moser 
(KAM) Q curves and generate a rich dynamics in quasi- 
integrable systems Q. In the context of soft walls, the 
sticky motion has been observed theoretically and exper- 
imentally in the one particle atom-optic billiard [l(| El 
and has shown to affect the quantum conductance in the 
soft wall microwave billiard [12| . 

Interacting many-particles systems in soft wall billiards 
is the next step to be studied. Such billiards are in- 
teresting not only from the fundamental point of view 
in nonlinear systems [l3| but also in many applications. 



*E-mail address: mbcims@fisica.ufpr.br 



Two recent examples can be mentioned: The effect of 
wall roughness in granular Poiseuille flow [Tij ] and how 
the confinement of the equilibrium hard-sphere fluid to 
restrictive one- and two-dimensional channels with soft 
interacting walls modifies its structure, dynamics, and 
entropy 0, [l6[ . It has been shown recently [ItJ that the 
origin of chaotic motion of two interacting particles in a 
onc-dimensional box is due double collisions which occur 
very close to the hard walls. These double collisions oc- 
cur when one particle is colliding with the ID wall and 
almost simultaneously collides with the other particle. As 
a consequence, the kind of motion generated close to the 
ID walls is essential for the whole dynamics inside the 
billiard. Therefore we expect that the softness of the ID 
walls will strongly affect the dynamics of the interacting 
particles. 

In this contribution we generalize previous results [l7l . 
[Lsj to the case of ID soft walls. The equivalence between 
the two interacting particles in the ID soft walls model 
with the motion of one particle inside a soft right triangu- 
lar billiard is shown. In this right triangular description 
the role of all important parameters from the problem be- 
comes clear. The interacting particles dynamics is stud- 
ied by varying the mass ratio 7 = mijm\ of the particles 
and the smoothness of the ID wall potential. The rea- 
son to use the mass ratio as a dynamical parameter is 
related to the generation of new materials in the field of 
nanotcchnology, where electrons may be confined inside 
a disk and can be affected by the surrounding material 
which composes the semiconductor [l9| . The composi- 
tion of the surrounding material changes the effective 
mass between particles 0, [2l[. We show here that a 
small "degree of softness (cr)" of the ID walls, strongly 
decreases the mean values of the finite-time Lyapunov ex- 
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ponents (FTLEs). The statistics of the distributions from 
the FTLEs has been studied in a number of physical sit- 



to Hamiltonian 
, and conserva- 



uations ranging from turbulent flows 22 
dynamics (in man y-p article system [22 
tive mappings [2~il. l25j]L In this work we use the FTLEs 
distribution over initial conditions, and the number of 
occurrences of the most probable FTLE, which has been 
proposed as an efficient quantity to detect small is- 
lands (dynamical traps) in phase space, to describe the 
qualitative and quantitative appearance of regular and 
sticky motion as a function of a . 

The paper is organized as follows. In Section [II] the 
model with soft ID walls used in this contribution is pre- 
sented. The description of this model in the right trian- 
gular is given in Section Mil Section IIVI shows a system- 
atic numerical study for the FTLEs, i. e., the number 
of occurrences of the most probable FTLEs as a func- 
tion of the mass ratios and of the softness parameter. 
Poincare Surfaces of Section (PSS) are used to show: a) 
The dynamics in the limit of hard ID walls and b) the 
rise of sticky trajectories. We end with the conclusions 
in Section fVl 



II. THE ID SOFT WALLS MODEL 

In this section we introduce the model used for the 
investigation of two interacting particles inside the ID 
billiard with soft walls. This model exhibit a continuous 
transition between soft and hard walls. The Hamiltonian 
considered is 



H = J2[T l + V l (q l )} + V mt =E, 



(1) 



where T,; = J^- (i = 1, 2) is the kinetic energy of particle 
Vint = Vo/r is the Coulomb repulsion between particles 
with r = \qi — q2\, and Vi{q{) is the potential energy from 
the ID soft walls, given by 



the walls looks very similar to the ID hard wall. This 
limit will be called here as the ID quasi-hard wall limit. 

The corresponding left and right force of the walls on 
particle i is 



Fj 



1 



1 



V2~7 



V2~7 



(3) 



The force of the walls has the Gaussian form, which in 
the limit a — ► approaches the 5-function. In this limit 
the corresponding potentials approach the hard 

walls. Rescaling the time by dr/dt = v / 2~E, the effective 
hamiltonian of the ID soft walls model is 



H=j=J2[Ti + Vi( qi ) 



V ir 



(4) 



where V int = V~o/r and Vi{q{) is given by Eq. ([2]) by using 
Fg instead Fq. The scaled potential intensities are Vq = 



H and F 



. These scaled intensities show the role of 



the total energy on the dynamics. Results of the present 
work, given for a combination of the scaled parameters, 
are valid for other energies keeping the scaled parameters 
constant. 
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t r / n Fa \ f (q l -d w \ ( q l + d w \\ 

VM) = y [erf (-^=-) - erf (— ^jj + F . 

(2) 

The first term on the right hand side of Eq. © represents 
the soft ID wall located at q = d w , while the second term 
represents the ID soft wall located at q = — d w . Here a is 
the parameter which quantifies the "softness" of the walls 
and Fq is the ID walls intensity. Fq is added such that 
the total energy is positive. Figure [T] shows the potential 
@ for different values of the softness parameter: a = 
5.0 x 10~ 3 (filled line), a = 5.0 x I0~ 2 (cross points), 
it = 9.0 x 10 -2 (dotted line). The ID soft walls are 
located at q w = ±1. As the softness parameter increases 
the walls become soft. An example of the ID soft wall is 
shown by the dotted line in Fig. [T] for a = 9.0 x 10~ 2 . As 
a approaches zero (see the filled line for a = 5.0 x I0~ 3 ) 



FIG. 1: Potential energy and corresponding forces for Fo = 
10.0 and a = 5.0 x 10" 3 (filled line), a = 5.0 x 10~ 2 (cross 
points), a = 9.0 x 10" 2 (dotted line). 

The model with smooth walls presented here differs, to 
our knowledge, from all models used in the literature. In 
our case, we assume the (5-function (in the limit a — > 0) 
for the force instead for the potential, as usual. The rea- 
son to do this is simple: when we use the J-function for 
the potential, then the corresponding force is not well 
defined. Our model makes it possible to study the classi- 
cal dynamics continuously in the transition from smooth 
to hard walls. Beside that, hard walls discontinuities in 
numerical simulations may introduce errors at each wall 
collision. 
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III. THE ID SOFT MODEL AND THE SOFT 
TRIANGULAR BILLIARD 

It is well know 0, [27], that the motion of three par- 
ticles on a frictionless ring with point-like interactions is 
equivalent to one particle moving freely inside the bil- 
liard, colliding elastically with the sides of the triangle. 
Let us start with the Hamiltonian of the three particles 

2 2 2 

on a ring, given by H B = >5 £ + & + W We as ~ 
sume that the elastic collisions take place at q\ = q%, 
92 — 93 and <j3 = qi + 1. These collision points define 
the hard-wall sides of the triangle. Using the orthogonal 
transformation [2(| (M = mi + 7772 + 7713): 
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92 



93 



for the three particles on a ring, the resulting Hamilto- 
nian is Hb = \i? + iy 2 + \z 2 . This is equivalent to the 
motion of one particle inside a triangular billiard with 
angles 



tana 



m 2 M 



tan/3 



miM 



tan 77 



m 3 M 



rn\m 3 V rn 2 m 3 y mim 2 

The point-like collision between particles 1 and 2 defines 
one side of the triangle at qi — q 2 = 0, and the collision 
of these particles with particle 3 defines the other two 
sides of the same triangle. For 7713 — > 00 (77 = 7r/2) we 
get the right triangular billiard which corresponds to the 
motion of two particles mi and m2 moving inside the ID 
box with hard walls. In this case the interaction between 
particles 1 and 2 is the point-like collision and the fixed 
particle 3 plays the role of the ID hard-wall, kept fixed at 
(73. These collisions with the ID fixed hard- wall can be 
represented by delta functions. However, the correspond- 
ing equations of motion are not well defined. Therefore, 
to describe analytically such problems we include in Hb 
the soft interactions between particles which, in a given 
limit, are expected to describe the point-like collisions 
and the ID box hard- walls. 

Firstly we assume that the interaction between par- 
ticles 1 and 2 is the Coulomb repulsion V12 = Vint = 
^o/| 9i — 92 1- In order to give an idea of the transi- 
tion to hard-walls in the triangle, we use the Yukawa 



potential V12 = V$e 



(-|9i-9a|/f) , 



|9i — 92 1- The interac- 



tion between particles 1 and 2 depends only on the rel- 
ative positions of both particles and v > is a param- 
eter which allows to change the interaction range from 
the smooth potential. For v — > the limit of short in- 
teractions (point-like collision) can be approached. Us- 
ing the orthogonal transformation (J5J) the above interac- 



tion is written as Vi 2 ( | y | ) = 



where 



fi 12 = mim2/(m.i-|-m2) is the reduced mass between par- 
ticles 1 and 2. Using V12 (|2/|), the one side of the triangle 
located at y = is now soft, as will be shown numerically 
later. 

Secondly, we assume that particles 1 and 2 interact 
smoothly with particle 3. In order to study the ID box 
case, the location of particle 3 is kept fixed at q 3 = ±d w 
and 7713 — ► 00. In this way particle 3 will play the role of 
the ID soft-walls. Here we consider the soft interactions 
Via(x, y) = Vi(qi) given by Eq. Due to the soft inter- 
action with particle 3, particles 1 and 2 can, in the ring 
description, interact with particle 3 on both sides, left 
((73 = — d w ) and right (q 3 = +d w ). In the limit 7773 — > 00 
the orthogonal transformation (f5j) is then reduced to the 
(x, y) plane 
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(6) 



V TO l m 2 



Rescaling the time by dr/dt = V2E and using the re- 
duced orthogonal transformation (|6|) we obtain, after 
straightforward calculation, the final scaled Hamiltonian 
Hb = Hb/(2E) = 1/2 in the right triangular description 
[using (i, y) -> (x, y)] 



H B = ^r + %r + V 12 {\y\) + ]T V l3 (x, y), 



(7) 



where 



Ui 2 (M) = V - 

\y\ 

V l3 (x,y) = F + 

\/2/ii2 

1m%o 




(-1) 



i+l, 



V m l m 2 



V2d v 
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erf 



V2/J12 

2m i g 



(-1) 



i+l. 



v / m 1 m 2 



2<T 



The total energy of the scaled problem is Hb = 1/2 and 
the effect of the real energy can be seen in the scaled 
potential intensities Vq = 



2E 



and F = |§. In 



the above Hamiltonian we observe the role of the cru- 
cial quantities of the problem, i. e. , masses mi, 777,2, the 
softness (v) of the interaction between particles 1 and 2, 
the softness (er) of the ID soft-walh, and the size d w of 
the ID soft walls billiard. 

Figure [5] shows the potential energy from ([7]) for two 
values of the smoothness a, and for two values of the 
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FIG. 2: (Color online) Scaled triangular potential energy 
(Vjjri) from the Hamiltonian with F = 10.0, Vo = 0.5 
and a = 5.0 x 10~ 3 in (a) and (c), and a = 2.0 x 10 _1 in (b) 
and (d). For (a) and (b) we have mi = 777,2 — 1, and for (c) 
and (d) we have mi = 3 and m,2 = 1. 



masses mi and 7712. We assume that V12 = ^b/M- First 
observation is that in the [x, y) plane the potential energy 
has a right triangle form, as expected. The left side of the 
triangle (y = 0) has a smooth form due to V12G2/I), which 
is the long Coulomb repulsion between particles 1 and 2. 
The other two sides of the triangle (with internal angle 



77 = 7r/2) can be soft or not, depending on a. These two 
sides are equivalent to the ID soft walls. For a = 5 x 10 -3 
and mi = m-x = 1 we nicely see in Fig.[^a) that these two 
sides approach to hard-walls, while for a = 0.2 these two 
sides are very soft, as can be seen in Fig. HJb). Figures 
[2jc)-(d) show the potential energy for the right triangle 
when the masses are changed (mi = 1 and m.2 = 3). 
Fig. djc) for cr = 5 x 1CT 3 and Fig. C^d) for a = 0.2. 
Above results show us that the motion of two particle 
inside the ID soft walls is equivalent to the motion of 
one particle inside the right triangle with soft walls. 

Some limiting situations can promptly be seen from 
Hamiltonian (0: (a) When F = (no ID walls) the 
two sides with internal angle 7/ — ir/2 of the triangle 
disappear and just the soft wall V12 ) remains. The 
Hamiltonian is separable and integrable since the x de- 
pendence in the potential energy disappears. This means 
that for two interacting particle (through relative coordi- 
nates) which are not bounded inside walls, the problem 
is regular, as expected; (b) The kind of coupling between 
coordinates x and y depends on the form of the inter- 
action used for V^x.y). Therefore, these two sides (or 
the ID soft-walls) determine if the dynamics inside the 
right triangle (or inside the ID box) is chaotic/regular; 
(c) When 777,2 — ► 00 (or mi — > 00) the x dependence in 
the argument of the error function (in Va{x, y)) can be 
neglected when compared to the y dependence, and the 
Hamiltonian ||7J) is separable again. This is also an ex- 
pected result, since when particle 2 (or 1) is too heavy, 
the motion should be regular again; (d) increasing the 
billiard size d w , the ID soft walls get apart. The same 
effect is obtained by decreasing the softness cr; (d) the last 
limit which we would like to mention, and the most in- 
teresting one, is the hard-wall billiard case, where a — > 
and v — ► (and Vq — > 00). In this limit the whole dy- 
namics can be explained in terms of rational/irrational 
values of a/7r [H, [29[. Using the Yukawa interaction 
in the ID hard walls billiard, we approached this limit 
numerically (y = 0.1) in a previous work [l7j . In [l7j 
it was also shown analytically the influence of the in- 
teraction between particles 1 and 2 to generate positive 
Lyapunov exponents. The softness of this interaction is 
essential to generate positive Lyapunov exponents, and 
it is possible to see that for the point-like collisions the 
probability to obtain positive Lyapunov exponents goes 
to zero. To study this limit analytically, the interaction 
potentials between all particles must be chosen appro- 
priately, and it is the subject of a future work. In this 
context it would be interesting to use the methodologies 
developed in [J, 0, 0, Q to study the limit a — > 0. 

We finish this section by mentioning that the Hamilto- 
nian ||7J) corresponds to the motion of one particle inside 
the right triangle suffering soft collisions at the walls. Al- 
though this is not, strictly speaking, a billiard motion as 
in the original sense (free particle inside a table), we will 
refer to it as the particle inside a soft right triangular 
billiard, due to the analogy shown in this section. 
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IV. RESULTS 

We investigate the non-linear behavior inside the ID 
soft billiard by determining the distribution P(A t ,j) of 
the finite-time largest Lyapunov exponents At as a func- 
tion of the mass ratio between particles 7 = 7712/mi. This 
investigation was done previously for the ID hard-wall 
case jirl . Il8j and the whole dynamics depends strongly on 
7. The FTLE is obtained by integrating two closed tra- 
jectories, computing the local LE after a time r = 0.1 and 
making an average over all the local LEs. The time r is 
chosen for the better convergence of the LE. We used the 
fourth-order Runge Kutta method with variable steps. 
The energy is conserved in all simulations by around 
10~ 6 . For quasi-integrablc system, the presence of broken 
KAM curves inside the chaotic sea in phase space, leads 
to trapping and 'sticky' trajectories [|[, affecting the con- 
vergence in the determination of the FTLEs, which de- 
pend now on the initial conditions. On the other hand, 
it implies that the distribution P(A t ,7), calculated over 
many initial conditions, contains information about the 
amount of regular motion (and trapped trajectories) in 
phase space [T3, HH HH HH HII • For ergodic system 
and for infinite times, the Lyap unov exponents do not 
depend on initial conditions [35j |. 



A. The ID quasi-hard limit (a = 5 x 10~ 3 ) 

Figure [3] shows (dashed line) the mean FTLE (A t ) as a 
function of the mass ratio 7 for a = 5 x 10~ 3 . This is the 
case of very small values of a (ID quasi-hard wall limit) 
and the ID soft walls look very similar to the ID hard- 
wall potential (sec Fig. [T|). The mean FTLE decreases 
from roughly ~ 1.15 to 0.54 in the whole mass ratio in- 
terval. This means that the dynamics is getting more and 
more regular, as expected, since 7 — * 00 constitutes an 
integrable limit with the heavy particle at rest. However, 
at 7 ~ 1.0 the mean FTLE increases, showing that the 
symmetry for equal masses increases the amount of irreg- 
ular motion. Figure [3] also shows (full line) results for the 
case of hard ID walls, published previously [Hj|. In that 
case, and therefore in the case observed here, the peak at 
7 = 1.0 is the consequence of the resonance observed in 
the limit of hard-point collisions between particles. The 
qualitative behavior of the FTLEs could be explained 
in [l8j with the help of some special periodic orbits from 
the Gauss map. For the purpose of this work we can ob- 
serve that both curves in Fig. [3J are in good agreement. 
This confirms that the model of ID soft walls presented 
here reproduces correctly, in the limit a — > 0, the ID hard 
wall case. However, small differences can be observed in 
Fig. [3J for 7 < 0.4,7 ~ 0.96 and 7 > 3.0. These differ- 
ences will be explained later, where we also show some 
Poincare Surfaces of Section (PSS) and discuss more de- 
tails of the whole dynamics for different values of the 
mass ratio. 

To analyze deeper the effect of the smooth potential 




FIG. 3: Mean values of the finite-time largest Lyapunov expo- 
nent calculated over 200 trajectories up to time t = 10 4 and at 
scaled energy E = 0.5, Vo = 1.0, F = 10.0, for a = 5 x 10~ 3 
(dashed line). This is the ID quasi- hard wall limit. For com- 
parison, the full line shows results obtained for the ID hard 
walls [131 . For each trajectory the largest FTLE is evaluated 
over 10 5 samples. 



on the particles dynamics, Fig. [?] shows the finite-time 
distribution of the largest Lyapunov exponent, P(A t ,7), 
for the (a) ID quasi-hard wall limit a = 5 x 10~ 3 and 
(b) the ID hard wall case from [Isj |. The gray points 
below the main curve are related to chaotic trajectories 
which were trapped for a while close to regular islands. 
Since both figures are quite similar, the ID quasi-hard 
limit also represents adequately the ID hard wall case 
concerning trapped trajectories. An interesting feature 
in Fig.0]is the change of the width of P(A t , 7) around the 
number of occurrences of the most probable A^ defined 
through [F^ ] 



0P(A*,7) 



dA, 



= 0. 



(8) 



For mass ratios close to 7 ~ 1.0, for example, many initial 
condition lead to different values of A t . In this region, 
A j has a minimum as a function of 7, which is a clear 
demonstration of the presence of 'sticky' trajectories. 

To determine the amount of 'sticky' and regular tra- 
jectories in phase space we follow A£ as a function of the 
mass ratio 7. In fact, we follow P(Aj), which is the nor- 
malized number of occurrences of the most probable LE, 
or the probability to obtain Af. This is shown in Fig. [5] 
When P(Af ) is large, a large fraction of initial condi- 
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FIG. 4: (Color online) Finite-time distribution of the largest 
Lyapunov exponent P(At, n /) calculated over 200 trajecto- 
ries up to time t = 10 4 , for (a) ID quasi-hard wall limit 
a = 5 x 10~ 3 , and (b) ID hard wall limit from [H]. With in- 
creasing P(At,7) the color changes from light to dark (white 
over yellow and blue to black). 



tions lead to the same A t and trapped trajectories are 
rare. For example, the maximum of P(A%) in Fig. 03 (a) 
[and (b)] close to 7 ~ 1.8, is the region in Fig.[4|a) where 
gray points below the main curve are rare. On the other 
hand, close to 7 ~ 0.96 we have a minimum in Fig. 03 (a) 
[and (b)], which is the consequence of the large dispersion 
around 7 ~ 0.96 in Fig. Ufa). Again the quasi-hard limit 
[Fig. 03 (a)] and hard wall [Fig. 03 (b)] agree very well. The 
fast variation of P(Af ) is due to statistical fluctuations 
in its determination over initial conditions. 

Figure 03 shows the PSS for 7 = 1.0, 1.8 and compares 
the ID quasi-hard limit [(a)-(b)] with the ID hard wall 
case [(c)- (d)]. Both cases are alike. The PSS is con- 
structed in the following way: each time particle 2 is lo- 
cated at the origin, and P2 > 0.0, then the point (qi,pi) 



FIG. 5: Normalized values of the most probable Lyapunov 
exponent Af for (a) ID quasi-hard case a = 5 x 10~ 3 and (b) 
ID hard wall case from uM- 



is recorded. Since this work is focused on the dynamics 
inside the realistic ID box problem, we have chosen to 
present results of the PSS in the ID box coordinates, in- 
stead in the right triangular coordinates (x,y). In fact, 
we have checked some PSS in the (x, y) coordinates and 
no additional relevant informations for the present work 
were obtained. It can be observed in Fig. 03 that for 
7 = 1.8 the PSS is almost filled with the chaotic trajec- 
tory, while for 7 = 1.0 some regular islands and forbid- 
den trajectories appear which induce the trapping tra- 
jectories (These trapped trajectories can not be seen in 
the resolution used in Fig. 03)- Therefore, in addition 
to the mean FTLE and the number of occurrences of the 
FTLE, also in the PSSs the quasi-hard limit and the hard 
ID wall case are almost identical, showing that our soft 
model correctly describes the hard ID wall limit a — ► 0. 
This will change drastically at next, when the softness 
parameter increases. 
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FIG. 6: Poincare surfaces of section for particle 1 and (a)-(b) 
ID quasi-hard limit and (c)-(d) ID hard wall case. Figures 
(a) and (c) for 7 = 1.0, and (b) and (d) for 7 = 1.8. 



B. 



The ID soft case a = 5 x 10" 



Finally we discuss the case of ID soft walls. The form 
of the ID wall potential is shown in Fig. [T] (see cross 
points for a = 5 x 10~ 2 ). Although the value of a is 
small, the effect on the particles dynamics is astonish- 
ing. Figure [7] shows the mean values of the finite-time 
largest Lyapunov exponents for a = 5 x (full line) 

compared with the quasi-hard ID wall limit (dashed line) 
from Fig. [3] The effect of the small soft potential is re- 
markable. The mean FTLEs decrease around 50% and 
the peak observed close to 7 = 1.0 (full line) is now a min- 
imum. A new minimum also appears close to 7 ~ 0.40. 

To understand better what happens, Fig. [5](a) shows 
the finite-time distribution of the largest Lyapunov ex- 
ponent for this case. Clearly we sec that, although the 
peak close to 7 = 1.0 still exist, the dispersion around 
the most probable FTLE is very large. This large disper- 
sion affects the mean Lyapunov exponents from Fig. [3] 
and a minimum occurs at 7 = 1.0. This strong disper- 
sion around 7 = 1.0 is also confirmed by the accentuated 
minimum of the number of occurrences of the most prob- 
able FTLE plotted in Fig.JSJb). A similar effect, but with 
smaller intensity occurs for 7 ~ 0.4. The mean FTLEs 
increase as 7 decreases, the amount of gray points below 
the main curve increases and the number of occurrences 
of the most probable FTLE also decreases. 

Many other values of the mass ratio 7 could be dis- 
cussed to analyze the appearance of island in phase space. 
Instead of doing so we would like to show some PSS in 
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FIG. 7: Mean values of the finite-time largest Lyapunov ex- 
ponent calculated over 400 trajectories up to time t = 10 4 and 
at scaled energy E = 0.5, Vq = 1.0, F = 10 for the ID soft 
wall with a — 5 x 10~ 2 (full line). For comparison, the full 
line shows results obtained for the quase-hard ID wall from 
Fig. El For each trajectory the largest FTLE is evaluated over 
10 J samples. 



order to discuss the dynamics of the interacting parti- 
cles in the ID soft walls. We start showing the case 
of 7 = 1.0, for which many trapped trajectories and 
islands are expected. Figure [Ha)-(c) shows the corre- 
sponding PSS. Since for the PSS we used the condition 
P2 > 0.0, particle 2 is always moving to the right at 
172 = 0. Because the interaction between particles is re- 
pulsive, we expect that particle 1 crosses the PSS more 
often for p\ < 0.0 than for pi > 0.0. This physical ef- 
fect is the origin of the asymmetry in Fig. [9ja) observed 
in p\. Additionally, due to P2 > 0.0, trajectories on the 
PSS with negative momentum (pi < 0.0) are less affected 
by particle 2 than trajectories with p\ > 0.0. Since the 
chaotic motion results from the simultaneous effect of 
the ID walls and the mutual interactions, we conclude 
that trajectories with p\ < 0.0 should be more regular. 
Therefore we expect more island for p\ < 0.0 than for 
pi > 0.0 [See Fig.^a)]. Observe that for p Y < 0.0 the is- 
lands extend themselves from q\ ~ —1.0 until q\ ~ —0.1, 
which is very close to particle 2, located at qi = 0. For 
Pi > the interval of islands goes from q\ ~ —1.0 until 
q 1 _o.4, very far from particle 2. Close to the left ID 
wall (<7i ~ —1.0), the momentum of particle 1 remains 
in the interval —1.0 < p\ < 1.0, while close to particle 
2 (qi ~ —0.1) we have chaotic trajectories with essen- 
tially only positive momenta (0.0 < p\ < 1.0) due the 
repulsion of particle 2. For —1.0 < p\ < 0.0 we have 
essentially regular trajectories since both particles are 
moving apart. The main island in Fig. [5{a) is related to 
the following regular trajectory: particle 1 is moving to 
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FIG. 8: (Color online) (a) Finite-time distribution of the 
largest Lyapunov exponent P(A t ,j) calculated over 400 tra- 
jectories up to time t = 10 4 , a = 5 x 10 -2 . With increasing 
P(At,7) the color changes from light to dark (white over yel- 
low and blue to black); (b) Corresponding normalized number 
of occurrences of the most probable Lyapunov exponent . 



the left with a small momentum which does not allow it 
to "penetrate" very much the ID soft wall. It "collides" 
with the ID soft wall and simultaneously is affected by 
the interaction force with magnitude comparable to the 
magnitude of the ID wall force. This is called here as the 
low-energy double collision, which in this example occurs 
between particle 1 (low-momentum), particle 2 and the 
ID soft wall. After colliding with the left ID wall, par- 
ticle 1 moves to the right with a small momentum and 
can not approach particle 2 very much due to the long 
Coulomb repulsion. Subsequently it changes its direction 
and goes back to the left ID wall. 

Another interesting property occurs very close to the 
"top of the ID soft wall potential" . We call the readers to 
attention that "top of the ID soft wall potential" in our 
one-dimensional billiard is the closest point that particles 
approach the turning point . In other words, it is the point 



where the "penetration" into the soft ID wall is maximal. 
When particle 1 has enough momentum (pi = —4.0) to 
"climb" ("penetrate") the left ID wall until the "top", it 
reaches the turning point on the left (around q\ <~ —1.1) 
with zero momenta [see Fig.[SJa)]. After that, particle 1 
returns to move to the right, accelerates (pi = +4.0) and 
travels until the other extremum at q\ ~ —0.2. This is 
a non-periodic chaotic trajectory with a behavior similar 
to the regular trajectory discussed above. The main dif- 
ference is that the regular trajectory has not enough mo- 
mentum to "climb" ( "penetrate" ) the ID wall potential 
until the "top" . In this example of the chaotic trajectory, 
the high-energy double collision occurs between particle 

1 (high- momentum), particle 2 and the soft ID wall. It 
is worth to mention that in the triangle description the 
above high- and low-energy double collisions correspond 
to the high- and low-energy particle collisions with the 
corners of the triangle which are located at y = 0. 

Figures [9^b)-(c) show a magnification of some regular 
regions from Fig.[5Ja), revealing the existence of the regu- 
lar islands. Figure[5]Jd) shows the PSS for 7 = 0.8, which 
compared to 7 ~ 1.0 presents a small number of islands. 
This is in agreement with results from Fig. [Hb), where 
the number of occurrences of the most probable FTLE 
has not a minimum for 7 ~ 0.8, and a small number of 
trapped trajectories is expected. 

We conclude this section saying that the softness of the 
ID wall decreases the mean FTLE. This effect is so strong 
that it is responsible for the small differences observed for 
7 ^ 0.4,7 ~ 0.96 and 7 > 3.0 between the quasi-hard- 
limit and the hard wall case [see Fig. ©]. For these 
values of 7 the number of occurrences of the FTLE from 
Fig. ([5]) has a minimum, and a larger amount of sticky 
trajectories is expected. This means that the softness of 
the ID wall has a stronger influence on space-phases with 
sticky trajectories. In addition we note that the above 
differences are stronger for 7 > 3.0 [sec Fig. This 
can be nicely explained using properties of the particle- 
particle frontal collision case. When a frontal collision 
of the particles occurs, then for 7 = m-z/mi > 3.0 the 
momentum of particle 2 is large enough and, after one 
collision, it continuous to move in the same direction as 
before the collision. Obviously this depends on the par- 
ticles energy, but in average, more frontal collisions are 
necessary to change the direction of movement of particle 

2 when 7 > 3.0. This property increases the amount of 
double collisions close to the ID wall, and consequently, 
the effect of the ID soft wall is more pronnouced when 
7 > 3.0, and the dynamics becomes more regular when 
compared to the ID hard wall case. 



C. Dependence on the softness parameter a 

In this section we analyze the mean FTLEs, (At), and 
the number of occurrences of the most probable FTLEs 
as a function of the parameter o. All simulations were 
realized over 400 initial conditions and for the mass ratios 
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FIG. 9: Poincare surfaces of section for a — 5 x 10 2 and for (a) 7 = 1.0, and their magnification (b) and (c). Figure (d) is for 
7 = 0.8. 



7 = 0.80, 1.0, 1.8,2.2. Figure [10] shows (A t ) for 10 differ- 
ent values of a in the interval 5 x 10" 3 < a < 5 x 10~ 2 . 
This is exactly the transition region between the ID 
quasi-hard limit and the ID soft wall. We observe that 
the main behaviour for all mass ratios is that the mean 
FTLEs decrease as a increases. This means that the de- 
gree of chaoticity decreases when the softness of the ID 
wall increases. 
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FIG. 10: Mean values of the finite-time largest Lyapunov ex- 
ponent calculated over 400 trajectories up to time t = 10 4 as 
a function of a, for 7 = 0.8, 1.0, 1.8, 2.2. 
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FIG. 11: Normalized number of occurrences of the most prob- 
able Lyapunov exponent A£ as a function of a. 



the dynamics is more ergodic-like for 7 = 0.8. This re- 
ally happens and can be confirmed visually comparing 
Figs. [Ha) and (d). Clearly the amount of regular island 
inside the chaotic sea is smaller for 7 = 0.8 [Fig. [9^d)] 
than for 7 = 1.0 [Fig.^a)]. 



V. CONCLUSIONS 



Figure QT] shows the number of occurrences of the 
most probable FTLE as a function of a for the values 
of 7 shown in Fig. QJJ For all mass ratios, but one 
(sec filled circles 7 = 1.0 ), P(Af) increases with a. 
This means that the dispersion around the most prob- 
able FTLE decreases when a increases, i. e., many initial 
conditions converge to the same Lyapunov exponent, and 
a more ergodic-like motion is expected when compared 
with other mass ratios. This is interesting seeing that 
for higher values of a the degree of chaoticity decreases 
but the dynamics becomes more ergodic-like. Compare 
for example the values of Af in Fig. [TTJfor a = 5 x 10~ 2 
and 7 = 0.8, 1.0. For 7 = 0.8 we have P(A?) - 0.2, and 
for 7 = 1.0 we have P(A£) ~ 0.1. So we expect that 



To conclude, we discuss the effect of physical realizable 
potentials (soft potentials) on the dynamics of interacting 
particles inside ID billiards. This contribution general- 
izes previous results [13, [H| to the case of ID soft walls. 
The ID soft walls are modelled by the error function with 
the softness parameter a. In the limit a — > the ID wall 
force on the particle is given by the ^-function. This 
allows us to study continuously the dynamics of inter- 
acting particles in the transition from soft to hard walls. 
The equivalence between the two interacting particles in 
the ID soft walls model with the motion of one parti- 
cle inside a soft right triangular billiard is shown. The 
role of all parameters from the model becomes clear in 
the right triangular description. Since the chaotic mo- 
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tion of interacting particles inside ID billiards is gener- 
ated by double collisions which occurs close to the ID 
walls [17| , we expect the influence of the soft ID walls on 
the particles dynamics to be very strong. Using the mean 
FTLE and the number of occurrences of the most prob- 
able FTLE, we analyze the dynamics of the interacting 
particles when the softness parameter changes from the 
ID quasi-hard limit (a = 5 x 10~ 3 ) to the ID soft wall 
case (cr = 5 x 10~ 2 ). We show that the ID quasi- hard 
wall limit agrees very well with results from the ID hard 
wall case, analyzed previously [l8| . When the softness 
parameter increases to a = 5 x 10 -2 , the mean FTLEs 
decreases substantiously (around 50%) when compared 
to the FTLEs from a = 5 x 10" 3 . Although both ID 
walls arc visually very similar, the dynamics of the in- 
teracting particles changes considerably. While regular 
islands and trapped trajectories are induced by the low- 
energy double collisions, the chaotic motion is produced 
by the high- energy double collisions. Double collisions 
are characterized by the simultaneous particle-particle- 
1D wall collisions. The rise of trapped trajectories is 
shown by using the number of occurrences of the most 



probable FTLE and the corresponding PSS. 

Results from the present paper strongly suggest that 
the transport of interacting particles and heat conduc- 
tion in physical devices [HI M, QJ, H, M, H, H, [H , 
which can be described by open billiard models, will sub- 
stantiously be affected by physically realizable wall po- 
tentials (soft walls). In such models, particles (heat) are 
injected at one open end of the billiard, and the efficiency 
of the transport depends how long (among other proper- 
ties) particles (heat) will need to reach the other open end 
of the billiard. Therefore, trapped trajectories induced 
by the soft walls and interacting particles, as shown in 
this paper, may substantiously increase the time spend 
by the particles (heat) inside the soft billiard, affecting 
the transport efficiency. 
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